Entanglement switching via the Kondo effect in triple quantum dots 
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We consider a triple quantum dot system in a triangular geometry with one ol the dots connected 
to metallic leads. Using Wilson's numerical renormalization group method, we investigate quantum 
entanglement and its relation to the thermodynamic and transport properties, in the regime where 
each ol the dots is singly occupied on average, but with non-negligible charge fluctuations. It is 
shown that even in the regime ol significant charge fluctuations the formation ol the Kondo singlets 
induces switching between separable and perfectly entangled states. The quantum phase transition 
between unentangled and entangled states is analyzed quantitatively and the corresponding phase 
diagram is explained by exactly solvable spin model. 



CO 



I 

C 

O 

o 



> 

in 
cn 

o 



PACS numbers: 73.63.Kv, 03.67.Mn, 72.15.Qm 



I. INTRODUCTION 

Creation of entangled states is essential for quantum 
computation and communication where most qubit op- 
erations cannot be performed through the manipulation 
of separable states Since the electron spin is a nat- 
ural two-level system the production and manipulation 
of spin qubit pair entanglement attracted much atten- 
tion 043 • Two electron spin qubits, each localized in 
one of two adjacent semiconductor quantum dots (QDs), 
can be coupled via the Hciscnberg exchange interaction 
due to virtual electron tunneling between the dots p|. 
The description of electrons by the spin degrees of free- 
dom alone is a simplification valid when the electrons 
are localized and the charge fluctuations arc negligible. 
In general, in any realistic solid-state device, spin en- 
tanglement is closely connected to the orbital degrees of 
freedom of the carriers which can be traced out if not 
measured @. Following the analysis of the use of en- 
tangled electron spin pairs in solid-state structures @, 
an intense activity was aimed toward understanding the 
physical mechanisms that produce spin-entangled elec- 
trons in mesoscopic conductors. The coherent manipula- 
tion of a single electron spin in a QD and the controlled 
correlation of two spins located in isolated dots have al- 
ready been demonstrated experimentally @. 

The interaction of qubit pairs with the environment 
is in general a complicated many-body process and its 
understanding is essential for experimental and theo- 
retical solid-state qubits |{|. Such a system is often 
represented by a two- level system, or spin- 1/2, inter- 
acting with an otherwise homogeneous, and often one- 
dimensional medium with gaplcss excitations. Some ver- 
sions of this model are equivalent to the Kondo model, 
motivating studies of ground state entanglement of an 
impurity spin with the conduction electrons and the role 
of the Kondo effect. This entanglement can be easily ex- 
pressed exactly in terms of the impurity magnetization 
H [T3| • It was found that the Kondo effect plays a signif- 



icant role in spin qubit double-quantum dot (DQD) (Io| . 




FIG. 1: Triple quantum dot system attached to the leads. 

Recently, the transport properties of triple quantum 
dots (TOP) in various configurations have been investi- 
gated [nj-yjl- In a serial TQD, different Kondo regimes 
can be sampled by measuring the conductance. The tran- 
sition from local-moment to molecular-orbital regime can 
be observed in the evolution of correlation functions [ll[ . 
In a system containing a central interacting dot C at- 
tached to the leads and two non-interacting dots A and 
B, there is a sizable splitting of the Kondo resonance 
[l2| and a two-stage screening of the magnetic moment 
was found. The inter-dot hopping introduces the Kondo- 
assisted transport and may induce a quantum phase tran- 
sition The possibility of a robust underscreened 
Kondo effect in TQD has also been recently discussed 

In this work we concentrate on thermodynamic and 
transport properties and their relation to entanglement 
for three electrons confined in three adjacent dots form- 
ing a ring; one of the dots is attached to metallic leads, 
Fig. [T] We study the case in which the dots arc singly 
occupied on average. The Kondo effect can switch on 
the entanglement due to the interplay between the inter- 
dot spin-spin correlations and various Kondo-like ground 
states. We find conditions for entangled and unentangled 
states between the dots. The results differ from those in 
the DQD configurations in which the Kondo interaction 
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reduces the entanglement between the DQD qubits (loj . 

The outline of the paper is as follows. Sec. II intro- 
duces the model and numerical method. In Sec. Ill nu- 
merical results are presented and summarized in a phase 
diagram and effective model in Sec. IV. 

II. MODEL AND METHOD 

We consider a solid-state qubit system built from a 
triangular TQD. We model it by the three- impurity An- 
derson Hamiltonian. The concurrence will be used as 
the measure of spin entanglement [ltj HH • We tune the 
parameters so that each dot is approximately singly occu- 
pied. The total spin is conserved in this system, thus the 
concurrence is related to spin-spin correlation functions. 

The Hamiltonian has three parts: H = Hd+Hi ea d+H t . 
The isolated dots are described by 

Hd = eid]di + U ^2 n %,<J n i,s + (1) 

i i 

+ E (*«4a- + h - c -)> ( 2 ) 

i<j.<J 

where = d\ a di a is the occupation of dot % € 

{A, B, C} with spin a. Only one energy level q = e 
is taken into account in each quantum dot - other levels 
are assumed to be much higher in energy. Parameter U 
describes the Coulomb interaction between two electrons 
in the same orbital. We use compact notation for the 
hopping integrals: t x = t A c, h = t AB , and t 3 = t B c, 
Fig. [TJ The dot C is coupled to the left (L) and the right 
(R) electrode, 

H t= ( V * c ika d c<r + h.c), ^ 

where V a is the tunnel matrix element between lead a 
and dot C. Electrons in the leads are described by a non- 
interacting Hamiltonian 

Hi ea d = ^2 (ekcl ka c akcr + h.c), (4) 

where with wavevector k is single-electron energy. d\ a 
(d^) and c^ fecr (c a ka) are the creation (annihilation) op- 
erators in the dot and in the lead, respectively. The leads 
have constant density of states p = 1/(2D), where D 
is the half-bandwidth of the conduction band. The hy- 
bridization strength is T = np (|Vl| 2 + |Vr| 2 ). Through- 
out this paper we assume Vl = Vr and t\ = fa . 

To properly describe the Kondo effect p7Hl9l |. im- 
purity problems need to be solved using nonperturba- 
tive methods, such as the numerical renormalization- 
group (NRG). The NRG [U Qj| consists of a loga- 
rithmic discretization of the continuum of states of the 
conduction-band electrons, followed by a mapping to a 
one-dimensional chain Hamiltonian with exponentially 



decreasing hopping constants. As a consequence of the 
logarithmic discretization, the hopping along the chain 
decreases exponentially, t n w A - ™/ 2 , where A is the dis- 
cretization parameter and n is the index of the site in the 
chain. This provides the opportunity to diagonalize the 
chain Hamiltonian iterativcly and to keep only the states 
with the lowest lying energy eigenvalues, since the energy 
scales are separated. Knowing the energy eigcnstates and 
eigenvalues, we can calculate thermodynamical and dy- 
namical quantities directly (e.g., spectral functions using 
their Lehmann representation). The method is reliable 
and rather accurate. Calculations in this work are per- 
formed with a discretization parameter A = 2, four val- 
ues of the twist parameter z, and the truncation energy 
cutoff of -E cu toff = IOwat, where ujm is the characteristic 
energy scale at the Nth NRG iteration. 

Since the NRG is a well-established numerical method 
with a relatively broad literature, we do not give the 
details of the implementation here but refer to Ref. 1171 
[20l - [22l | . In this work we use the NRG Ljubljana code [18| . 



III. RESULTS 

A. Quantum entanglement 

Our goal is to analyse the entanglement of a qubit sub- 
system attached to external reservoir. A similar study of 
DQDs has shown that at low temperatures and low mag- 
netic field, the Kondo effect acts as a source of entangle- 
ment destruction 0, [13] ■ On the contrary, in the present 
triple quantun dot system, the Kondo effect enables the 
entanglement, as shown in the following. 

Quantum entanglement of two-qubits in a pure state 
\AB) = ottl tt) + a n \ U) + c^fl it) + a u \ U) can 
customarily be quantified by the Wootters concurrence 

M 

Cab = 2|atta:^ - a n a^\. (5) 

Two qubits arc fully entangled, Cab = 1, if they are in 
one of the Bell states: ± lit) or |tf) ±\U)- 

In TQD, the qubits are coupled to a fermionic bath, 
therefore the two-qubit system AB is in a mixed state 
and appropriate generalization for concurrence is given 
by the Wootters formula [l(|. Moreover, in the system 
considered there are also charge fluctuations, thus the 
state may be outside the simple manifold of spin degrees 
of freedom. Nevertheless, the concurrence can still be 
consistently defined for such system. In the presence of 
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axial spin symmetry, it can be given as @, [l(]| 
max(0,C ti ,C||) 



Cab = 



Pn + P\\ 



(6) 



a 



n 



2\{StSs)\-2j(P\Pl)(PiPi\ (7) 



C\\ = 2\(SlSi)\-2j(PiPi)(PiPl), (8) 



P 



n 



P\ 



pt pi I pi pt 

A B ' A B I i 



pt pt i pi pi 

A B r A r B 



(9) 
(10) 



where 5 4 + = (S*~)^ = is the electron spin raising op- 

erator for dot i = A, B and P? = rij CT (1 — ?i i(T ) is the pro- 
jection operator onto the subspace where dot i is singly 
occupied by one electron with the spin a. P|| and Pj-j, are 
probabilities for the spins to be aligned in the same (par- 
allel) and opposite (antiparallel) directions, respectively. 
This formalism is also applicable for the concurrence be- 
tween other dots, e.g. Cqa (concurrence between the dot 
C with the dots A) and Ccb (concurrence between the 
dot C with the dots B). The expectation values can be 
easily obtained from the NRG. 



B. Thermodynamic quantities and transport 

In order to better understand the nature of the cre- 
ation of the entanglement studied here we additionally 
calculate: 

(1) The temperature-dependent impurity contribution 
to the impurity magnetic susceptibility |17l | 



Ximp(P) 



k B T 



(11) 



where S z is the z component of the total spin of the whole 
system while (...) means the thermodynamic expectation 
values. The first expectation value refers to the system 
with dots, while the second (with the subscript 0) refers 
to the system without dots; hb is the Bohr magneton, g is 
the g-factor, and k B is the Boltzmann constant. It should 
be noted that the combination k B Tximp{T) / \gn B ) can 
be considered as an effective moment of the impurities. 
(2) The impurity contribution to the entropy (l7T | 



<S'imp(P) — 



(E -F) (E- F) 



T 



T 



(12) 



where E = (H) = Tr[Pexp(-P/fc B T)] and F = 
-k B T\nTr[exp(-H/k B T)}. 

(3) Thermodynamic expectation values of various op- 
erators such as the on-site occupancy (rii), local charge 
fluctuations ((Jn^) 2 ) = (nf) — (n^) 2 , and spin-spin cor- 
relations (S, ■ S,-). 



(4) The electronic transport between the leads through 
dot C, making use of the Meir-Wingreen formula [23| 



G(T) = GqttT 



duj 



A c (u,T), (13) 



where Go = 2e 2 /h is the conductance quantum, / is the 
Fermi function, and Ac(u>,T) is the spectral function on 
the impurity C. 



C. Numerical analysis 

The analysis of spin-spin correlation functions demon- 
strates the singlct-triplct transition between the dots A 
and B. The two regimes are separated by a first order 
quantum phase transition (level crossing) as a function 
of t2/t\ which determines the ground state of the iso- 
lated TQD (l2T - fl4| . The same quantum phase transition 
has been previously found in a double quantum dot sys- 
tem modeled as a pure spin system |24j |. Entanglement 
properties of a qubit pair formed by spins on adjacent 
quantum dots are closely related to the spin-spin cor- 
relations @ therefore transitions between different spin 
configurations play a crucial role. 

We first describe the strongly correlated regime with 
U = 10r at temperature below the Kondo scale (essen- 
tially we are in the T = limit). Each of the dots is 
almost perfectly singly occupied. It should be noted that 
contrary to some double QD systems here exact single oc- 
cupancy of identical dots in general cannot be achieved 
for a common value of e since the system is asymmetric 
and, furthermore, it is not particle-hole symmetric due to 
the lack of bipartiteness of the lattice. Expectation val- 
ues (...) in the concurrence formula Eq. ([6]) correspond to 
the thermal equilibrium of the system and consequently 
(«!?^ Sg') =0. In vanishing magnetic field, the concur- 
rence formula (JS|) simplifies further to 

Cab = max {0, -2(S A • S B )/(P n + P\\) - 1/2} . (14) 

Therefore, the concurrence is significant with increased 
spin-spin correlations in the range —3/4 < (Sa ■ Sb) < 
-1/4. 

In Fig. [2](a) we show the concurrence between the ad- 
jacent dots. As expected, there are two regimes of inter- 
est, depending on the inter-dot couplings, similar to the 
result of a recent analysis of thermal entanglement of iso- 
lated TQDS [H . On one hand, for t 2 > t 2c we get perfect 
entanglement between the dots A and B, which are an- 
tifcrromagnctically coupled, Fig. [^b). Critical intcrdot 
coupling t2 C is equal to t\ for decoupled quantum dots, 
r = 0, and is slightly renormalized due to the coupling to 
the leads [H, O, H3 ■ F° r increasingly large coupling t 2 , 
the spin-spin correlations diminish due to charge fluctu- 
ations, Fig. [UJd), which reduce the probability for single 
occupation of the dots, Fig. H£c). However, the concur- 
rence remains constant, Cab = 1- This means that two 
electrons would form a perfectly entangled qubit pair, if 
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io" 3 10" 2 10"' 10° 

t 2 /D 



FIG. 2: Concurrence, spin-spin correlations, probabilities for 
parallel and antiparallel spin configuration, charge fluctua- 
tions, entropy Si mp /k B , susceptibility k B Txim P {T) / (g^ B ) 2 , 
and conductance at temperature T = I0 _6 D, as a function of 
t 2 /D for U/D = 0.1, Y/D = 0.01, e = -U/2, h/D = 0.01. 



extracted one from the dot A and the other from B. The 
dot C is in the Kondo regime with electrons in the leads 
and it is completely decoupled from the other dots. Sim- 
ilar behavior, known as a dark-spin state, occurs in TQD 
ring in presence of an in-plane electric field due to decou- 
pling of the spin in one of the dots [|| . On the other hand, 
for t% < ti, the concurrence is zero and ferromagnetic 
correlations between dots A and B dominate, forming 




io" 3 io" 2 io" 1 io 1 

T/D 



FIG. 3: Concurrence, spin-spin correlations, probabilities, 
and charge fluctuations at zero temperature as a function of 
T/D for fo/ti = 0.5; other parameters as in Fig. [2] 



an effective 5=1 impurity at low temperatures which 
undergoes partial Kondo screening, yielding a residual 
uncompensated spin-I/2. 

In addition to the concurrence, spin and charge cor- 
relation functions, abrupt changes also occur in thermo- 
dynamic properties. In Fig. Hfe) we show the entropy 
change from S- lmp /kB = In 2 to zero and the suscepti- 
bility from Ximp = 1/4 to 0, Fig. [^f). Finally, at the 
transition the linear conductance changes from zero to 
G = Gq. This is expected due to the decoupling of the 
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dot C from the entangled pair AB - the conductance is 
unity as in the case of a single Anderson impurity. 

The transition at t 2c ~ ti is explained by the level 
crossing of the corresponding doublet eigenstates of iso- 
lated TQD [| |, GJ, El- However, for t 2 much lower 
than t\ full-entanglement between the dots A and B can 
be restored by a sufficient increase of the effective Kondo 
coupling of the dot C to the leads, leaving the dots A 
and B in the singlet state. In Fig.[3|a), we show how the 
concurrence Cab abruptly changes from zero to unity 
with increasing dot-lead hybridisation T. This is an al- 
ternative way of inducing the same quantum phase tran- 
sition, as evident from the entropy, susceptibility and the 
conductance which all exhibit behavior analogous to the 
results presented in Fig. [2] The main difference is in the 
fact that dots A and B remain consistently decoupled 
from the rest of the system for the whole range of T, as 
is seen from the spin-spin correlations signaling triplet or 
singlet state for T < T c and T > T c , respectively, while 
{Sa ■ Sc) ~ due to the increase of Suq with increasing 
T/D. Charge fluctuations 5n A remain low, as expected. 



IV. DISCUSSION AND SUMMARY 

Results presented in Fig. [2] and Fig. [3] are typical ex- 
amples. In order to test the robustness of the results with 
respect to charge fluctuations on the dots, we performed 
a detailed numerical analysis summarized in the phase di- 
agram in the parameter space (U/T, t 2 /t\), see Fig.[4ja). 
The moment of the dot C, (S c ), is represented as a con- 
tour plot (dashed lines). It is clear that for larger T (i.e., 
smaller U/T) the moment is diminished which partially 
contributes to the renormalization of t 2c . 

The main reason for the renormalization can, in fact, 
be understood from an effective spin model by tracing out 
the empty and doubly occupied states, and representing 
the conduction band using a single spin S d 

M, 



represen ting 

BIlSliI 



H c ff = Ji(Sa + Sb) ■ Sc + J2&A ■ Ss + Jd&c ■ &d, (15) 

with the exchange couplings Ji 2 = Tt\ 2 /U . First we 
consider a simple TQD system with Jd = 0, i.e., with- 
out spin Sjj. The eigenstates are two doublets and a 
quadruplet (which always lies the highest in energy). 
The ground states for the z-component of the total spin 
S z = 1/2 are given by 



\Di) = 



I Da 



(I UIb) - I Ufa)) ® I tc 

V2 



(16) 



2J tttg) ® ||c) - (| UIb) + I Uts)) g I tc) 
V6 



The first state \D\) is formed from the singlet state be- 
tween the dots A and B, for which the spin of the dot C 
is decoupled from the other dot spins 
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FIG. 4: (a) Phase diagram in the (U/T,t2/ti) plane for 
ti/D = 0.01 and fixed U/D = 0.1. Full red line (connect- 
ing calculated tac/ti, bullets) separates Cab = 1 (the Kondo 
phase) and Cab = regions between two ground state con- 
figurations. Dashed lines represent constant values of local 
moment (Sq). (b) Phase diagram in the (Ji/Jd, yJT/Ji) 
plane for an effective pure spin model, as an analogue of (a) 
with full line representing yJ^c/Ji- 



The second state |D 2 ) is constructed from the triplet 
states between the dots A and B with S z = 1 and 0, 
separately. The spin-spin correlation functions between 
the dots A and B for the above two states are 



(Di|S A .S B |Di) = -3/4, 
(D 2 \S A ■ S B \D 2 ) = 1/4. 



(18) 



Therefore, for the states \D±) and \D 2 ) the spin-spin cor- 
relations between the dots A and B are anti- and ferro- 
magnetic, respectively. The eigenvalues are 



E 2 



-12%/U, 
-4(4^ 



4)/u. 



(19) 



(£>i|S c ■ S A |Di) = (£>i|S c • S B |Di) = 0. 



(17) 



From the energy separation E\ — E 2 = J\ — J 2 it is clear 
that J 2c — J\ is the critical value separating entangled 
from uncntanglcd states which explains results for large 
U/T. 

Phase diagram presented in Fig. HJa) can qualitatively 
be understood also in the regime of lower U /T where 
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t2c is significantly lower that t±. A simple analysis of 
the effective Hamiltonian (|15jl for Jd > 0, i.e., for TQD 
coupled to an additional singly occupied quantum dot D, 
reveals a crossing of levels representing different ground 
states where the dots A and B form triplet or singlet 
configurations, if the interaction Jd is below or above 
some Jdc, respectively. Critical value Jdc can be given 
analytically, 

. (Ji- J 2 )(Ji + 2J 2 ) 

= ^j" • ( 2U ) 

In Fig.U^b) we show the phase diagram (Jx/Jd, \fJ2JJ1) 
where the full red line represents the separation between 
entangled and uncntangled regimes. The topology of the 
dots is shown as an inset. Here Jd/ J\ plays the role of an 
effective Kondo coupling of the TQD to the leads, Jd 
T/U. The dots A and B for small J2/J1 but large Jd/ Ji 
thus exhibit perfect entanglement when dots C and D are 
forced to couple into a singlet, which can be considered as 
an example of the entanglement monogamy concept [25| . 
However, due to the substantial renormalization of the 
local moment (S^) in the limit of small values, U/T < 4, 
the application of the pure spin model is not fully justified 
there. This is the reason for different behavior of critical 
t2c/t\ oc U/T and \J J% c / J\ cx \J J±/Jd in this regime. 

In summary, the spin entanglement of electron pairs 
in a triangular TQD, one dot being attached to non- 
interacting leads, is quantitatively analyzed in the regime 
of two competing many-body effects, namely the Kondo 
effect and the direct exchange interaction. In contrast 
to DQD configurations in which the Kondo effect sup- 
presses the entanglement [Toj], in this case, the Kondo 



effect - in which a single spin in the dot C is screened 
by conduction electrons in an attached metallic lead - in- 
duces entanglement between the spins in the dots A and 
B. There are two regimes of interest, depending on the 
ratio of the interdot couplings t%/t\. On the one hand, 
for t2 ^ ti, spins in the dots A and B are always max- 
imally entangled and the central dot C is in the Kondo 
regime with conduction leads, thus the entanglement is 
switched on due to the Kondo effect. For stronger in- 
teraction regime U/T ^> 1 this result is expected as a 
consequence of the quantum phase transition observed 
in TQD with triangular topology [HI, [l3| . Our analysis 
reveals that the switching between unentangled and per- 
fectly entangled qubit pairs A and B persists also in the 
regime of substantial charge fluctuations, either 5n\ or 
Sn,Q. In the regime of low U/T there is a strong renor- 
malization of t2 C <C t±, which can be understood within 
a simple analytically solvable pure spin model where the 
TQD is anti-ferromagnetically coupled to an additional 
singly occupied quantum dot. The induction of the sin- 
glet formation C-D acts as switching of the entanglement 
between the dots A and B in a similar manner as the for- 
mation of the Kondo singlet in the original TQD system 
analyzed numerically. 
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